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Abstract 

Principal components analysis coupled with non-parametric 
bootstrapping is introduced with an example as a powerful 
tool to help visualise, analyse and decide on comparative or 
single non-deterniinistic LCI/LCIA results. Decision support is 
provided by adding non-parametric bootstrapping (NPB) to 
the GAIA plane, which is a special case of principle comments 
analysis (PCA) built around the Promethee multicriteria deci¬ 
sion aid model. 

In addition to an easy to analyse visual presentation of otherwise 
complex and cluttered numerical results, PCA with NPB is equally 
at ease will all data uncertainty formats used to date in LCA. 
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1 Introduction 

The analysis of LCA results can be made difficult in practi¬ 
cal situations by the quantity of data, the multiple units types, 
the various media to which substances outflow, judgmental 
values to be applied, etc. Analysis is also made difficult by 
the uncertainty on background and foreground data. As a 
direct consequence of this, LCA results can be difficult to 
communicate to non specialists. 

In the "classical" LCA model this situation is handled by 
not taking uncertainty into account (possibly getting some 
insight afterwards through sensitivity analysis, for instance) 
and by producing a single index for each alternative. Doing 
so, comparing alternative multidimensional data sets pro¬ 
duced by alternative LCI and LCLA becomes straightforward. 
This position, however, is now evolving in many ways. Not 
only is blind confidence in the "classical" LCI and LCLA 
models and their results decreasing, but also uncertainty in 
the input data is more often reckoned with. Thus LCA is 
steadily evolving into a non deterministic model. 


Depending on the practitioner, non-deterministic LCAs will 
not always be performed up to the index calculation step. If 
they are, however, results are no longer multidimensional 
but only uncertain. Common statistical techniques can be 
used to document the index uncertainty and assess the rel¬ 
evance of the comparison of two indexes (Steen). If an in¬ 
dex is not calculated, the practitioner is faced with multidi¬ 
mensional and uncertain impact indicators or aggregated 
flows sets, thereby missing ways to get a synthetic view apart 
from a tedious analysis of tables like Table 1. 

Extracting synthetic information from multidimensional (or 
"multivariate") data sets is the purpose of a family of data 
analysis techniques. Some of them have already been used 
in LCA. Multivariate factor analysis is mentioned by (Huele) 
and principal components analysis (PCA) has been used by 
(Le TEno, 1997) on deterministic LCA results. In both works 
mentioned, these techniques are used to visualise "the whole 
thing" and analyse its inner structure. Also, (Mareschal, 
1988) presents a special case of PCA which can be used to 
provide direct support for multicriteria decision making. 

The purpose of this paper is to present example results of 
viewing, analysis and decision support with non-determin¬ 
istic LCA output data through the application of PCA cou¬ 
pled with non-parametric bootstrap (NPB) as an uncertainty 
analysis technique. Examples have been produced with 
MatLab, a commercial numerical computation package. A 
brief mathematical background is given in the appendix. 

2 Viewing Multidimensional and Uncertain Data 
Sets with PCA 

2.1 PCA with certain data 

For illustrative purposes, we will use comparative LCI/LCIA 
results from (Gemis) (—» Table V). Most criteria are precise 
values except for the last four ones which are intervals of quali- 


' Table 1 is a mixture of environmental and economical criteria, as will most 
likely be the case in the final analysis of LCA results. 
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Table 1: Alternatives x Criteria table resulting from multiple LCI-LCIA 
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4 
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13 
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0.29 
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0 
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0 
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0 
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tative values. As for now, only their mean value will be consid¬ 
ered. We will deal with their uncertainty in the next paragraph. 

The purpose of PCA is to produce a planar view of the ta¬ 
ble, so that most of the information is preserved. To do so, 
each alternative (line) A. is considered as a point whose co¬ 
ordinates are given by the value on each criterion (column) 
(—> Figure 1). A plane on which the projection of the A.'s 
will preserve most of the initial shape of the cloud of A.'s is 
first identified. Then, correlation coefficients between crite¬ 
ria are projected as segments starting from the plane's origin. 
The computed PCA plane for Table 1 is shown on Figure 2. 



Fig. 1: Alternatives A, (matrix lines) as a point in the criteria’ space 
(matrix columns) 


Alternatives are drawn as stars and criteria as segments start¬ 
ing from the origin (0,0). Horizontal and vertical scales are 
not relevant. 

Significance of the vertical and horizontal axis is in terms of 
synthetic criteria driving the comparison. It can be derived 
from the length of the criteria's projections on each axis. 
Thus, alternatives are ranked according to their direct con- 
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Fig. 2: PCA plane of mean of Table 1 
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tribution to the greenhouse effect on the horizontal axis (cri¬ 
teria C,, C, and C n ) and the more we move left, the better 
the alternatives perform. There appears to be a cluster of 
rather similar alternatives in that direction (A n , A p , A |3 , 
A |4 , A |7 ). Logically enough, alternatives A 4 and A s (thermal 
electricity) are positioned at the extreme right of this axis. 
At the same time, the more right, the better the alternatives 
perform on criteria C,,, C [3 and C )4 , related to soil pollution 
and use. Obviously trade-offs will have to be found between 
the two groups of criteria, most probably among alterna¬ 
tives belonging to cluster A,, A,, A 3 , A s , A 1S and A 16 , which 
also perform moderately well on other gas emissions (C„ 
C,, C 5 , C 7 ), or alternatively to the A m , A |S , A |9 , A ? and A, 
group, if these gas emissions are not that important to the 
decider. Also, we find that C s , as a short segment, is not a 
very discriminating criterion. This can be verified in Table 1, 
where C s is not very different between alternatives. Concern¬ 
ing costs (C 9 and C ]0 ), they appear to be correlated to green¬ 
house gas emissions. Since we also know that C 10 (external 
cost) is calculated from CO, emissions in (Gemis), this should 
be no surprise. Hence, C 10 is somehow redundant with C 4 or 
C h . Internal costs (C 9 ) are correlated less frankly. They are, 
however, directly opposed to hazards (C p ). Again, some level 
of compromise will have to be found in a final analysis. 

We see that a rather straightforward but nevertheless deep 
understanding of the initial data set (—>• Table 1 ) can be gained 
from simple geometrical considerations. The same kind of 
analysis can be performed on a single LCI/LCIA, for instance 
with processes as alternatives and flows as criteria. This is 
partly what is done in (Huklk), as well. For a detailed de¬ 
scription of PCA and the foundations to PCA result inter¬ 
pretation, readers are referred to (Murtacjh), which also 
provides C source code or (Lebart). 

Now, how will our conclusions hold if data uncertainty of 
Table 1 is taken into account ? 


2.2 PCA with uncertain data 

In case data are uncertain A/s are no longer single points 
but families of points contained in (hyper)cubes in the crite¬ 
ria coordinates system. This is illustrated in Figure 3. Prob- 



Fig. 3: Alternative A ; with interval data on all three criteria 


ability distributions, fuzzy numbers or intervals apply along 
the (hyper)cube's edges and indicate the known distributions 
of values. Crisp numbers also apply. To deal with such data, 
PCA may produce zones of confidence around each projec¬ 
tion of A using either analytical techniques or simulation. 

A very common approach is to proceed by simulation using 
a re-sampling technique called non-parametric bootstrapping 
(NPB) (Lebart). With NPB, a deterministic PCA plane is first 
computed using mean values. Then, values from within the 
(hyper)cube around each alternative are sampled according 
to their distributions and projected on the PCA plane, pro¬ 
ducing "clouds", the hulls of which indicate how variable 
alternatives positions are. 

To proceed, Table 1 is turned into Table 2. When data are 
intervals (last four criteria), only lower and upper bounds 
are sampled with a 0.5 probability each, as no information 
are known on values within the bounds. The first 10 criteria 
in Table 1 are precise values. What we do for the purpose of 
this article is to add uncertainty according to a normal dis¬ 
tribution centred on the mean value and extending +/-10% 
on each side. Adding uncertainty is a parametric boot¬ 
strapping that should not be used when a non deterministic 
LCA is performed since empirical distributions or intervals 
on the data are known and can be used instead of more or 
less arbitrary ones. 200 additional matrix are calculated and 
plotted on the initial PCA plane computed on mean data. 
Approximately 30 additional matrix is a minimum (Lebart). 

PCA of Table 1 is shown in Figure 4. As can be seen, alter¬ 
natives are more or less spread over the plane - A |s and A w 
being very condensed and A 10 very widespread. As of the 
clusters of alternatives identified in the deterministic PCA, 
A n , A p and A,, are too overlapped to be considered as sig¬ 
nificantly different; A |7 is a little bit more separate and much 
more condensed. As for A |4 , it can be considered as separate 
from the other alternatives in the cluster. 



Fig. 4: PCA plane of Table 2 with bootstrap 
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3 Supporting Decision Making through PCA 

As illustrated on the examples above, a usual PCA visually 
provides relevant information to assist in the decision mak¬ 
ing process: 

• Synthetic criteria can be identified through projection of the 
criteria segments on the vertical and horizontal axis, thus 
giving a better understanding of the global factors at play 

• Highly correlated criteria can be identified and possibly 
simplified into fewer criteria 

• Highly influent criteria and negligible ones can be iden¬ 
tified through their length on the PCA plane 

• Clusters of similar actions can be located and possibly 
eliminated from the analysis 

Additionally, using the non parametric bootstrapping tech¬ 
nique, more certainty can be gained as: 

• The (non)-overlapping of two alternatives' hulls shows 
how much sense their ranking makes; alternatives which 
do not separate enough can be eliminated 

• All available data are visualised, not only mean values 

• A largely spread alternative which may deserve investment 
into a more precise data collection is clearly identified 

In a decision context, alternatives are plotted with a view to 
identifying the best performing one(s), and not only to locate 
them as function of their input and output flows, as previ¬ 
ously. As discussed in the Figure 2 comments, final decision 
depends on a compromise between criteria, but this compro¬ 
mise is not explicitly nor fully handled by PCA. This is where 
multicriteria decision aid (MCDA) models have a role to play. 

The Promethee model (Brans) is the one used in the EQuity 
model for an LCA of construction products (Le Teno, 1996) 
as an interpretation help tool. For each criterion, Promethee 
is asking the user to provide: 

• A weight 

• Whether a minimum or a maximum is sought 

• A preference function, indicating how the preference 
varies as function of the difference between alternatives 
on the given criterion 

Six types of preference functions are available, allowing lev¬ 
els to be set as well as linear or sloppy growth in preference. 
Example preference function ("type 6" in Promethee) is 
shown on Figure 5. 

Promethee includes a variant of PCA called " GAIA" as a visual 
decision support tool. Decision relevant information is re¬ 
quested from the user in a clear and user-friendly manner. Prior 
to performing deterministic PCA as usual, the decision matrix 
is preconditioned in a way that integrates these preference infor- 



Fig. 5: Promethee example preference sygmoid function ("type 
6 ", parameter is o) applied to the GWP indicator, a and b are 
alternatives. 


mation except for the weights (—> Appendix). In addition to 
classical alternative and criteria plotting, a "decision vec¬ 
tor" is drawn pointing in the direction of the "best alterna¬ 
tives overall" (as function of the user's criteria weights). As 
far as our knowledge, except for the decision vector orienta¬ 
tion, uncertainty in the data is not handled in GAIA. What 
we present here is an example generalising the use of 
bootstrapping to deal with uncertainty in GAIA, hence pro¬ 
viding a more complete decision support than simple NPB 
PCA or deterministic GAIA. 

The NPB GAIA plane generated with data from Table 2 and 
preference information given in Table 3 is shown in Figure 6. 
As can be seen with alternatives A lg and A„ in Figure 4 or A 10 
in Figure 6, clouds are sometimes made out of disconnected 
"stacks” of dots, and not continuous spreads. This is more a 
(slight) visual annoyance than a real problem. In all cases, 
hulls are the only relevant concept, regardless of the distribu¬ 
tions of initial data within their bounds being continuous or 
not. The positions of the alternatives have changed a little from 
what they were in Figure 4, due to preference information. 
Given other preference information, it should always be ex¬ 
pected that the picture will be altered to precisely reflect this 
information. As for the criteria, their relative "lengths" have 
changed since, contrary to normal PCA, preference functions 
vary from one to the other. In the case drawn, previous com¬ 
ments still hold. Supplementary information is available 
through the decision vector: Alternatives A2, A3, A8 are close 
to the best compromise direction (according the user's prefer¬ 
ences) and could be chosen for further investigation or public 
debate, for instance. 

4 Some Comments on the Modelling of Uncertainty 
in LCA 

On a general basis, we think that resorting to crisp (mean, 
etc.) values in LCA/LCIA, when irrelevant, is a dangerous 
convenience since it may lead to unstable conclusions, 
hence ruining the current LCAs and LCA's predictive cred- 
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Tabic 3: Preference information on Table 2 criteria in the Promethee format 


Criteria 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

Weight 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Min / max 

min 

min 

min 

min 

min 

min 

min 

min 

min 

min 

min 

min 

min 

min 

Preference 
function type 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

Parameter (o) 

0.36 

0.32 

0.12 

466.18 

2.56 

2.97 

1.46 

0.03 

4.38 

2.90 

1.36 

0.49 

1.38 

1.34 


1 . 5 ! 



— 2 *- 1 - 1 - 1 - 1 - 1 - 1 - 1 - : - 1 - 1 

- 2.5 -2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 2.5 

Fig. 6 : GAIA plane of Table 2 with bootstrap 

ibility and public reconnaissance at the same time. Hence, 
providing strong support for uncertainties in LCA is a 
priority. 

As of the present situation, depending on the mathemati¬ 
cal format chosen for uncertainty, alternative data will be 
intervals (Le T£no, 1996; Heijungs), fuzzy numbers or 
probability distributions (Coulon, Steen). It is worth men¬ 
tioning, however, that uncertainty in data is also a multi¬ 
dimensional concept that may have various aspects (Klir), 
including: 

• Fuzziness-data are vague 

• Non-specificity-information are lacking 

• Dissonance-conflicting information is available 

Also, as shown in Table 4, all data uncertainty mathemati¬ 
cal formats do not reckon with all uncertainty aspects. Hence, 
even though various non-deterministic formats are suggested 
for LCA, they should not be considered as mutually exclu¬ 
sive, but rather as complementary. Ideally, non-determinis¬ 
tic LCA should be able to cope with all three types of uncer¬ 
tainty mentioned, hence with all uncertainty formats used 
so far, either directly or through a common meta-format 
including them all. 


Tabic 4: Uncertainty types modeled by some uncertainty models 


Uncertainty type model 

Fuzziness 

Non¬ 

specificity 

Dissonance 

Intervals 

1 

Yes 


Fuzzy sets 

Yes 

Yes 


Possibility distributions 


Yes 


Probability distributions 



Yes 


5 Conclusions 

Principal components analysis (PCA) coupled with non-para- 
metric bootstrapping (NPB) is introduced with an example 
as a powerful tool to help visualise, analyse and decide on 
comparative or single non-deterministic LCI/LCIA results. 
Decision support is provided by adding NPB to the GAIA 
plane, which is a special case of PCA built around the 
Promethee multicriteria decision aid model. In all cases, 
though PCA with NPB is equally at ease will any mixture of 
uncertainty formats which, aside its graphical nature, is a 
big advantage of this technique. 

It is hoped that viewing, analysis and decision support tech¬ 
niques as the one presented here will foster the turning of 
LCA into a non-deterministic model. 
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Appendix 

Mathematical foundations and necessary demonstrations are 
omitted for the benefit of a simpler and almost " ready to code" 
description of PCA and GAIA. Readers are referred to (Lebart), 
(Murtagh) and (Mareschal, 1998) respectively for details. 


Principal components analysis 

Let R be an alternatives x criteria (n,k) matrix. Columns of 
R are centred around their average (centering) and reduced 
to unit standard deviation. Let X be the resulting matrix. 
Let , the first two highest eigenvalues of matrix X'X (in 
decreasing order) and u,, their corresponding eigenvectors. 
(u p u 2 ) defines the first PCA plane on which projections of 
alternatives and criteria will be made. Coordinates (A' and 
C) of alternatives A. and criteria C respectively in the (u p u,) 
plane are obtained as follows: 

A] =(X.u i ,X.u 1 ) 

C] =(x'x.u t /J^,x , x.u 1 /Jx^) 


GAIA PCA plane 

Let R be the decision matrix. R is transformed into matrix 
<t> with the following equation: 

<!>(/,7) = - l -E { p j < A i' h ) - P i<*■ A , )} ’ 

with A. = alternative number i, b an alternative belonging to 
A, the set of all alternatives A, and P.(.,.) the preference func¬ 
tion associated to criterion j. Then, proceed with d> as with 
matrix X above. To get the decision vector, simply plot a 
segment from the origin to the projection of vector e defined 
below on the PCA plane. 



w. is the weight of criterion j and ej the null vector with 
component j set to 1. 
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